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Abstract 


The first steps in developing a methodology for spacecraft control-structure interaction 
(CSI) optimization are identification and classification of anticipated missions, and the develop- 
ment of tractable mathematical models in each mission class. A mathematical model of a generic 
large flexible space platform (LFSP) with multiple, independently pointed rigid payloads (repre- 
sentative of “Class 11” CSI missions) is considered. The objective here is not to develop a general 
purpose numerical simulation, but rather to develop an analytically tractable mathematical model 
of such composite systems. The equations of motion for a single payload case are derived, and 
are linearized about zero steady-state. The resulting model is then extended to include multiple 
rigid payloads, yielding the desired analytical form. The mathematical models developed clearly 
show the internal inertial/elastic couplings, and are therefore suitable for analytical and numerical 
studies. A simple decentralized control law is proposed for fine pointing the payloads and LFSP 
attitude control, and simulation results are presented for an example problem. The decentralized 
controller is shown to be adequate for the example problem chosen, but does not, in general, 
guarantee stability. A centralized dissipative controller is then proposed, requiring a symmetric 
form of the composite system equations. Such a controller guarantees robust closed-loop stability 
despite unmodeled elastic dynamics and parameter uncertainties. 


Outline 


• CSI Mission Classification 

• Mathematical Model of LFSP/ Single Payload 

• Extension to Multi-Payload Systems 

• Decentralized Control Law 

• Simulation Results for Space Station/Articulated Payloads 

• Centralized Dissipative Controller 

- guarantees closed-loop stability 

- requires symmetric form of eqns. of motion 

• Concluding Remarks 

182 


H juttlii mum hi mu i 



CSI Mission Classification 


Anticipated missions involving control-structure interaction (CSI) can be broadly divided 
into four classes. Class I missions are those which require fine-pointing of the overall spacecraft 
as well as vibration suppression. There are no articulated substructures, and because of small 
elastic and rigid-body motion, the modeling and control problem is essentially linear. There may, 
however, be actuator and sensor nonlinearities. Examples of this mission class include large space 
antenna concepts such as the hoop-column and wrap-rib antennas. The controller must satisfy 
the performance specifications, which include fine-pointing (rigid plus flexible rotational motion), 
reflector surface distortion, and defocus errors caused by feed and base elastic motion. The control 
system must be capable of performing in the presence of significant elastic motion. 


• Class I : Single-Body Flexible Spacecraft (Linear) 

- fine-pointing and vibration suppression 

- surface shape distortion and defocus errors 



hoop-column antenna concept 


wrap-rib antenna concept 


183 


CSI Mission Classification (cont.) 


Class II missions represent an extension of Class I missions with articulated payloads 
mounted on the flexible spacecraft. This important class of missions includes many types of 
low-earth and geosynchronous platforms. The composite system consists of a large flexible space 
platform (LFSP) to which a number of rigid or elastic appendages are mounted. The objectives of 
the controller are precision attitude control of the spacecraft, fine-pointing of each of its payloads, 
and vibration suppression. The problem is still linear, but the elastic and rigid-body motions 
between bodies are coupled. 


• Class II : Flexible Spacecraft with Articulated Appendages (Linear) 

- includes rigid and clastic appendages (payloads) 

- fine-pointing of all rigid and elastic components 

- vibration suppression to improve payload performance 
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CSI Mission Classification (cont.) 


Class III and IV missions are essentially the nonlinear counterparts to the Class I and II 
missions previously discussed. Class III missions will require large angle manuevers of the entire 
structure (without articulated appendages) while simultaneously minimizing the effects of elastic 
motion induced by nonlinear internal couplings and external disturbances. Although this problem 
has been studied numerically for specific flexible spacecraft, generic models do not exist. 

Class IV missions have the inherent difficulty of the nonlinear Class III problems, but in- 
clude the additional complexity of the articulated appendages of their linear Class II counterparts. 
Such missions may require large- angle maneuvering of the flexible structure while simultaneously 
and independently pointing various payloads to their own respective targets in the presence of 
elastic motion. There is rigid/elastic coupling between all of the various components. Further- 
more, payloads may be repositioned, or the spacecraft reconfigured by commanded motions of 
the various appendages. Examples of this mission class include LFSPs with articulated payloads 
and robotic arms allowing translational degrees of freedom. 

The purpose of this paper is to develop a generic mathematical model and control laws 
for an LFSP with multiple articulated payloads representative of Class II missions suited to CSI 
optimization. 


• Class III : Single-Body Flexible Spacecraft (Nonlinear) 

- Precision large angle maneuvers for retargeting or tracking 

- Vibration suppression during maneuvers 

• Class IV : Flexible Spacecraft with Articulated Appendages (Nonlinear) 

- Large angle rotational motion of spacecraft and appendages 

- Robotic manipulators allowing spacecraft reconfiguration 
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Mathematical Model for LFSP/Single Payload 


In developing the equations of motion, the single payload composite system will first 
be described in terms of its various coordinate systems, position vectors, etc. The equations of 
motion for each component comprising the composite system, including LFSP rigid and elastic- 
motion as well as payload rigid-body motion, will then be presented. By removing the constraint, 
force between the LFSP and payload, these equations will then be coupled, and presented in 
matrix form. This generic model will then be extended to the multi-payload case. 

A large flexible central body with a single attached articulated rigid payload is consid- 
ered first. The mass/inertia properties of both bodies are assumed known. The modal solution 
(frequencies and displacement vectors) for n modes of the flexible body (LFSP) is also assumed 
known. This modal solution is usually obtained using finite element codes such as NASTRAN or 
EAL, wherein the effects of the payloads are taken into account by representing them as point 
masses and inertias on the LFSP model. The present analysis requires the values of LFSP mode 
shapes and mode slopes at the points of application of forces and torques, respectively. Both are 
required at the payload attachment point. The attitude control system for the LFSP is assumed 
to utilize reaction control system (RCS) thrusters, and control moment gyro (CMG) torques. The 
payload pointing mechanism is idealized as a three-axis point torque acting at the attachment 
gimbal. The objective is to control the payload line-of-sight and LFSP attitude in the presence 
of elastic motion. 

Overall Configuration 


• LFSP with single rigid payload 

• Control objectives: 

- attitude control of LFSP 

- fine-pointing of payload 


• Data required for model 

- payload and LFSP masses and inertias 

- modal frequencies 

- mode shapes for forces (e.g.: RCS thrusters) 

- mode slopes for moments (e.g.: CMG moments) 
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Coordinate Systems 


The LFSP has two coordinate systems associated with it. One, designated the o-system, 
is fixed to the nominal center of mass (CM), and in conjunction with the nominally zero Euler 
angles (<£ 0 , 9 0 , ?/> 0 ), describes the rigid-body orientation. The other coordinate system, designated 
the e-system, is fixed to the payload attach point (“point P”) on the LFSP, thus rotating and 
translating with the rigid and elastic LFSP motion at that point. This system is nominally parallel 
to the o-system axes. The orientation of the elastic axis system (e-system) with respect to the 
o-system axes is described by Euler angles (<; f> e , # e , ip e ) (i.e., the rotational elastic displacements), 
which are nominally zero. The origin of the LFSP o-system is located in the reference inertial 
coordinate system (denoted the n-system), by the position vector R no = Xh\ + Y112 + Zh 3 . The 
payload attachment point P for the undeformed spacecraft is located from the LFSP CM, the 
o-system origin, by position vector R op = x p 61 +V P 62 + 2pO S . 

The rigid payload is described by a body-fixed p-system originating from the attachment 
point, but fixed to the payload. The e- and p-system origins are coincident. Orientation with 
respect to the elastic axes is described by Euler angles (<^ p , 9 p , </> p ). These are the payload pointing 
mechanism “gimbal angles”. The desired nominal orientation is defined by user-specified gimbal 
angles. Thus the angular orientation of the payload may be determined from the orientation of 
the p-system with respect to the e-system; the e-system with respect to the o-system (due to 
elastic deformation); and the orientation in inertial space of the o-system. Expressing a vector 
in terms of components in the various coordinate systems requires transformations using a 1 - 2-3 
Euler sequence of rotations. The angular velocity vectors (body rates) of the various coordinate 
frames are u; 0 , u? e , and u? p for the LFSP rigid frame, LFSP elastic frame, and the payload frame, 
respectively. The body rates must be transformed into Euler rates prior to linearization. 

• inertial reference: n-system 

• LFSP (rigid) body axes: o-system 
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Nonlinear Model 


The nonlinear equations of motion for the composite system are developed next. The 
bodies are each assumed to have constant mass and mass distribution, and to possess only rota- 
tional freedom with respect to each other. All terms in the equations represent inertial quantities 
which are to be expressed in the o-system, nominally parallel to the inertial frame. It should 
be noted that w ^ ” denotes differentiation with respect to the inertial frame, and so angular 
velocity cross-products appear, whereas an overdot “ * ” indicates a local time derivative inside a 
particular frame. Appropriate coordinate transformations in terms of Euler angles, using a 1-2-3 
sequence, are required. 

First consider the LFSP equations of motion. Beginning with the translational equation 
of motion, writing Newton’s law gives equation (1). Here m 0 is the LFSP mass, F p is the internal 
constraint (reaction) force exerted on the LFSP by the payload, and there are n p external forces 
F 0 i acting on the body (e.g., RCS thrusters). 

Next the nonlinear rotational equation of motion of the LFSP is given by equation (2). 
I 0 is the LFSP inertia matrix about the CM, T p is the pointing system reaction torque, and 
there are moments M c ; acting on the LFSP (e.g., CMG moments). The position vector 
R op = x p 6\ + y p o 2 + z p ds locates the payload attachment point, and the position vector r 0 ; = 
x 0 i 6 i + y 0 i02 + ^0x63 locates the i th external force F 0 i acting on the LFSP. Here the “prime” 
notation with R4 p an d r oi indicates the (3x3) matrix cross-product form. 

LFSP Eqns. of Motion 

• LFSP translation: 

** F 

m oR-no = ^ p ^ \ ^oi 

• LFSP rotation: 

n,\t np 

Io^’o H" ^ lo^o “ Hop x Fp — Tp “f- ^ " M 0 i ~h ^ ^ r Q { x F 0 j (2) 

i=l t = l 
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Nonlinear Model (cont.) 


The LFSP elastic motion is assumed to be decoupled from the LFSP rigid-body motion 
(valid for small motion), so that inertial/elastic coupling terms arise from payload interactions 
with the platform at the attachment point. Thus the elastic motion of the LFSP may be de- 
scribed by equation (3), the modal equation of motion. The (n x 1) vector q contains the modal 
coordinates, where n is the number of modes used to describe the elastic motion. The modal 
coordinates and their corresponding time derivatives will comprise 2 n components of the system 
state vector. The matrix D is a diagonal matrix with terms 2 where and are the damp- 
ing ratio and frequency, respectively, for elastic mode i. Matrix A is a diagonal matrix of lo \ . 
The $ s consist of mode slopes at the point where a moment is applied, or mode shapes at force 
application points. The modal matrices are dimensioned (n x 3). The right-hand side of the 
expression represents modal forces and moments. 


LFSP Elastic Motion 


ti M Tip 

q 4- Dq -f Aq = -3> u F p - #oT p + ^ ^ 3>uoiF 0 i (3) 

i=i ;=i 
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Nonlinear Model (cont.) 

Next consider the payload equations. The payload CM location, point P*, for the de- 
formed spacecraft with respect to the inertial reference is given by vector R np * in (4) below, The 
acceleration of the payload CM is sought. Consider time derivatives of the position vector R np * 
with respect to the inertial frame, that is, take of the right-hand side of R np * • This gives 
the velocity of the payload CM, R np *. Taking another time derivative, again including terms 
due to the relative rotations between coordinate systems, yields the inertial acceleration of the 
payload CM, R np » • The only external force acting on the payload is the constraint force at the 
attachment point. Writing Newton’s law gives the nonlinear translational equation of motion, 
equation (5). 

Lastly, the nonlinear payload rotational equation of motion (about the payload CM)_is 
given by equation (6). I p is the payload centroidal inertia matrix for the p*-system, and M p is 
the total moment about the payload CM. Taking appropriate time derivatives and characterizing 
the moment M p acting on the payload results in equation (7). 


Nonlinear Payload Equations of Motion 

• Payload CM translation: 


R-np* — R-no + Rop 4" 4 Tp* 


m ptt„p* — Fp 


• Payload rotation: 

x d , 

Ip + Wp 


(4) 

(5) 
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Linearized Model 


Consider now transformation of all the inertial quantities to o-system components, lin- 
earization of the nonlinear equations, and the introduction of perturbations (to nominal values). 
Note that the o-system is nominally parallel to the inertial frame. Upon linearization, system 
variables may be thought of as being comprised of a nominal value plus a perturbation quantity, 
for example some variable y = yo+6y. The prefix “6” will be used to denote departure values from 
the equilibrium (nominal) state. The nominal attitude of the LFSP is a Q = (cf) 0 ^9 0y V>o) T = 0, the 
nominal angular elastic deformation at the payload attach point is = (<£ c , $ e ?^e) T — — 

0, and the nominal gimbal angles are a p = (<f> p ,O p ,ip p ) T = (<^ p , 6 p , ^ P ) T . Furthermore, the modal 
coordinates used to describe the elastic motion are nominally zero, as are the forces and moments 
acting on the system. 

The angular velocities of the various coordinate systems may be expressed in general as 
u) — uj 0 4 - Su), However, the nominal angular rates are zero, and <5u; = <$a, with a = (<f> , 9 , t/>) t , 
the Euler angle vector. The elastic deformation angular rate may be written in terms of the mode 
slopes and modal coordinates as follows: <5a e = 


• “Variable” = “Nominal” + “Perturbation”, y = yo + Sy 

• Nominal Values: 

- LFSP attitude, a Q = (</> 0 , 0 o , rfo) T = 0 

- modal coordinates, q = 0 

- gimbal angles, a p = (</> p ,0 p ,^ p ) T 
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Linearized Model (cont.) 

The linearized perturbation equations of motion are given below. The prime notat ion — 
denotes the matrix cross-product form. C ep is the transformation from the e-system to the 
system for the nominal gimbal angles a p = (<^> p , 0 p , */> p ) T . Also, I po = Cj p I p C ep is the payloacF 
centroidal inertia matrix I p (expressed in the p*-system) transformed into the o-system. I p |p 
premultiplies the payload inertial angular acceleration. 


• LFSP translation: 

n F 

m 0 5 R no = -5F p 4- ^2<5F 0 i (S) 

!=1 

• LFSP rotation: 

" Af 71 F 

1 0H0 = -KpWp - <^Tp + J2 6M °' + r oi 6F oi * (9) 

7=1 1=1 

• LFSP elastic motion: 

n m n F 

<5q + D<5q + A5q = -$ u T p - $<^T P + $o 0 i<5M 0 i + $uoi<5F oi (10) 

7=1 i=l 

• Payload CM translation: 

m p {«R„o- [Rip + tCjpl-p.)' ]«a 0 + [*3'-(Cj p r p .)'#J']<q ^ 

-(Cj p r p .)'Cj p <a p } = «F p 

• Payload rotation : 

Ipo(«o + + Cj p <5a p ) = <5T p - (Cj p r p .)'«5F p (12) 
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Linearized Model (cont.) 

The final equations of motion can be coupled by substituting for the constraint force 
d)Fp given by equation (11) into the remaining equations (8, 9, 10, and 12), thus eliminating the 
expression for payload translation. The LFSP translational, rotational, and elastic equations, 
and the payload rotational equation may be expressed in the form of equation (13). Here A is a 
(4 x 4) block matrix of overall dimension (n + 9) square. Matrices B and C are essentially null 
except for n diagonal terms corresponding to the elastic motion. The vector Sr} is (n + 9) x 1, 
and is given by equation (14). 

The expression for matrix A is given by (15). Various terms comprising matrix A are 
given by equation (16). The expressions for matrices B and C are given by equation (17). The 
matrices D and A were given previously. The mathematical structure of this generic model 
includes the internal couplings between payload rigid-body motion and LFSP rigid and elastic 
motion. 


• Model form: 


A8t j + BSi) -f C 5t) = Etui 


• Vector of system variables: 


Stj — (5(A,l,^,<^o, 0 q 


n i 0p i i ) 


(13) 


(14) 


• Matrices A, B, and C: 


-(?77 0 + m p )I 3 x3 

— 7T7 p R/ 

m p 4> 

-m p r'Cj p 

^^pH'op 

Io - ?77pR4 p R' 

m p R; p 4> 

— m P Ro p r'Cj p 

rrip&u 

-m p $ u k' 

Inxn 4" 

-m p # u r'Cj p 

m p f' 

I P o - m p r'R' 

Ipo^J + m p r'<& 

(I P o - m pf' r')Cj p 


where: 

f' = (Cj„r p .)' R' = (R; P + F') $ = (*! - T'*J) ( 1C ) 



06X6 


0gx6 

B = 

DnXn 

C = 

■A-n Xn 


0 3 x3 _ 


0 3 x3 _ 


193 


Linearized Model (cont.) 

Now consider the right-hand side of equation (13). The vector of control inputs is given 
by (18). Matrix E is made up of np block-columns corresponding to forces <5F 0 ;, n m block- 
columns for moments (5M 0 ;, and a single block-column multiplying the payload torque. The 
overall dimensions of matrix E are (n + 9) X (3 np + 3 n^f + 3), and E as a partitioned matrix is 
given by equation (19). Each of the three partitioned matrices comprising E are given in (20). It 
is clear that E x corresponds to external force inputs to the LFSP, E 2 to external moment inputs, 
and finally E 3 corresponds to the payload gimbal torque. 


• Vector of control inputs: 


6u = (6F^,...,6Fj-6M 


r T 

*ol * ’ 


; <5Tp) T 


• Partitioned form of matrix E: 


(IS) 


E - [ Ei( n+9 ) x(3nF ) , E 2 ( n +9)x(3n M ) > E 3 ( n+9)x3 


(19) 



’ 13x3 

^3 x 3 

r 


r ol 

r 7 

onp 


Ex = 

^uol 

• • * uonp 

, E 2 = 


-0 3 x3 

03x3 - 



03x3 

O 3 X 3 


O 3 x 3 

I 3 x3 

Ljx3 

, e 3 = 

_ 13x3 

^Ool 

' ■ ^(?onM 


-$o 

03x3 

0 3x 3 - 


- I 3 X 3 - 


( 20 ) 
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Decentralized Control Law 


The matrix equation (13) may be put into state variable form by allowing x x = 6r], as 
given in (14), and x 2 = Sr/. Thus one_obtains the expressions in (21), which may in turn be 
written as equation (22). The matrices A and B are given by equation (23), with each sub-block 
of A and B dimensioned (9 + n) square. 

The approach considered herein is a decentralized control strategy using simple propor- 
tional plus rate feedback with collocated actuators and sensors for both payload pointing and 
LFSP attitude hold. Only applied moments are considered, with LFSP attitude control effected 
by moment input <5M 0 , and payload pointing accomplished by an applied girnbal torque <5T p . 
The sensed LFSP attitude and attitude rate at the CMG location are given by (24). The LFSP 
attitude is the sum of the central body rigid angles, and the elastic rotations at the CMG location. 
The sensed attitude and attitude rate of the payload are given by (25). The payload attitude is 
the sum of the central body rigid angles, the elastic rotation angles at the payload attach point, 
and the girnbal angles. 


• State-space Form 

Xl = x 2 

x 2 = A -1 [Eu — Bx 2 — Cxi] 
x = Ax + Bu 



0 I 


0 

A = 


B = 



_ — A -1 C -A -1 B 


1 

> 

i 

H 
1 


( 21 ) 

( 22 ) 

(23) 


• Sensor Outputs 
— LFSP attitude and attitude rate: 
6a oa = <5ao + $£,<5q 

(24) 

6ao g = 5k 0 + $ J Q Sq 


— Payload attitude and attitude rate: 

<5a pg = Sa 0 + + Cj/a p 

(25) 

<5a pB = 6a 0 + + C ep <5a p 
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Decentralized Control Law (cont.) 

The decentralized control laws for LFSP attitude control and payload pointing are given 
by equations (26) and (27), respectively. Consider the LFSP attitude control system. If the 
gain matrices K ao and K ao are chosen as in (28), the resulting closed-loop system (assuming no 
coupling) would have a damping ratio of p Q , and a bandwidth of rad/sec. Similar expressions 
for the payload gimbal controller may be written by replacing A 22 with A 44 as in equation (29), 
giving a closed- loop system with damping ratio p p and bandwidth 


• Control Law 

- LFSP attitude: 

<5M 0 = — [K ao <5a os + K^a os ] (26) 

- Payload Gimbal Torquer: 

<5T P = -[K ap <5a ps + K* p <5a ps ] (27) 

- LFSP Gain Matrices: 

K a 0 = -A 22 diag(f2“) ^ 

K Ao = -A 22 diag(2/> 0 Q 0 ) 

- Payload Gain Matrices: 

K a p = — A 44 diag(fip) ( 29 ) 

K& p = -A 44 diag(2pp^p) 
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Extension to Multi-Payload Case 

The equations of motion for multi-payload configurations may be written in the form of 
equation (30). The vector 6rj is dimensioned (6 + n + 3n p ) x 1, and is given by (31). The vector 
of force and moment inputs <5u is dimensioned (3n^ + 3 Rm + 3 n p ) x 1, and is written as in (32), 
where 6F 0 ;, <5M 0 i, and <5Tj are external LFSP forces and moments, and payload gimbal torques, 
respectively. The coefficient matrices appearing in equation (30) are given in Appendix A. 


• matrix form of eqns.: 


A <5r? + B Sri 4- CSrj = E6u 
• vector of system variables: 


(30) 


Sr] — <5( A , V, Z ; 0 O , $o, V-’o i 5 i i ••• ■> Q n 5 » V’l > ••• > 4 > n p i @n p i ^Pn p ) (31) 


• vector of control inputs: 


<5u = 6( Fjj, ... , Fj nr ; M£, 


M 1 • T 1 

? M on M ’ A 1 ’ 


T \T 


T 1 ) 

’ x n D / 


(32) 
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Simulation Results for an Example Problem 

In order to demonstrate the modeling and control methods developed, the dual-keel space 
station ISS04 reference configuration with individually pointed payloads was considered. Mass, 
inertia, and other parameters for the station are summarized in the table below. Many different 
simulation scenarios were considered in exercising the space station model, however only a two- 
payload study comparing the effects of end-mounting versus CM-mounting will be presented 
herein. Attitude control moment input to the station was provided by a three-axis CMG near 
the station CM, while three-axis gimbal torquers acted at each of the two payload attachment 
points. In this study, the “worst-case” station rigid-body initial conditions were chosen based on 
the most severe anticipated maneuver, re-boosting to higher orbit. Initial modal coordinates were 
obtained from the steady-state elastic deformation resulting from the firing of RCS thrusters. The 
decentralized controller must drive the station rigid-body angles and payload pointing errors to 
zero in the presence of elastic motion. 


ISS04 Configuration j 

j 
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Simulation Results (cont.) 


The 12 vibrational modes that were actually used in the simulation studies, and their 
corresponding modal frequencies, are tabulated below. These 12 modes were chosen because 
they represent major truss-structure bending and torsional deformations, as opposed to mere 
appendage modes. The elastic modes were assumed to have a damping ratio of 0.005. Payload 
data are also summarized below. Given are the payload attachment points, and the masses and 
inertias. The first rigid payload was mounted on the lower keel, and the second payload on the 
upper keel. For the control systems, the desired closed-loop damping ratios were assumed to be 
0.707, and the desired bandwidths were 10.0 rad/sec and 0.1 rad/sec for the two payloads and 
the space station, respectively. 


Model Data 


mode # 

frequency 

(rad/sec) 

7 

1.37236 

8 

1.39753 

9 

1.50451 

20 

2.01425 

21 

2.24925 

22 

2.68481 

36 

4.09370 

37 

4.74973 

38 

4.87077 

39 

5.86502 

40 

6.51971 

41 

7.28311 


Modes Included in Composite System Model 


Payload Data 



payload 1 

payload 2 

1SS04 gimbal node 

3005 

3002 

position vectors 

x = 8.3333 

x = 9.1667 

from CM (ft) 

y = 5.4167 

V = -60.202 


* = 236.223 

* = -142.708 

payload mass 

100.0 

100.0 

(slug) 



payload inertias 



(slug-ft 2 ) 




1250.0 

1250.0 

t'V 

0.0 

0.0 

/*. 

0.0 

0.0 

Iyy 

1250.0 

1250.0 

lyi 

0.0 

0.0 


1500.0 

1500.0 
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Simulation Results (cont.) 


The steady-state (limit cycle) station rigid-body pitch and pitch rate encountered in a 
reboost maneuver simulation study were 0.02187 rad (1.25 deg) and 0.0008727 rad/sec (0.05 
deg/sec), respectively. The present simulation used these values as initial angular displacements 
and velocities about each of the station axes. The initial modal coordinates were calculated by 
assuming the space station had reached steady-state conditions after the the firing of the four 
RCS thrusters in the x - direction at 75 Ibf each. This assumed maneuver serves only to provide a 
reasonable initial defo rme d state. There are no thruster inputs after the start of the simulation. 
The payload gimbal Euler angles and the remaining states are initially zero. The desired gimbal 
angles for both payloads are ai i2 = 0.1(1, 1, 1) rad. It was desired to return the station to its 
nominal orientation, i.e. the Euler angles (<f> 0 , d 0 , t/’o) = 0, while driving the payload pointing 
error to zero, in the presence of nonzero initial conditions and any induced elastic deformation. 
Two cases were compared as a result of studying the equations of motion: an end-mounted case 
with nonzero distance between the gimbals and the CM, r p * — (0, 0, 6) T ft, and a CM-mounted 
case where r p * = (0, 0, 0) T . 


• Initial Conditions 

- each axis: angle = 1.25 deg, angular rate = 0.05 deg/sec 

- initial elastic deformation due to thruster firing (steady-state) 

- initial gimbal angles are all zero 


• desired gimbal angles: a p i = a P 2 = 0.1 (1, 1, 1) rad 

• station stabilization, simultaneous payload pointing 

• two cases: 

- end-mounted payload 

- CM-mounted payload 
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Simulation Results (cont.) 

The results given below are for Payload 2 on the upper keel, about the pitch axis (the 
x and 2 motions are analogous, as are the results for the payload on the lower keel). The figure 
shows the “base” motion, that is, the sum of the space station rigid-body orientation angles 
and the elastic rotation at the payload attachment point, for an end-mounted payload. This 
motion must be compensated for by the payload gimbal torquer so that the payload remains on 
target. The large rigid-body motion about the y-axis masks the much smaller elastic motion. 
The absolute payload pointing error for the end-mounted case is also shown below. The absolute 
pointing error is the sum of the base motion and the payload gimbal angles, and should be driven 
to zero by the payload controller. 

The inset figure compares the absolute pointing error for the CM- and end-mounted cases. 
Note that the absolute error for the end-mounted case is significantly more oscillatory than the 
CM-mounted case. The end-mounted case has a maximum overshoot of -136.49 arcsec. This 
response remains oscillatory with a magnitude on the order of 10 arcsec. The CM-mounted case 
yields the desired error response (i.e., perfect second-order system response) about the y-axis 
(and x and z axes as well) because the payload rotational motion is completely uncoupled from 
the LFSP motion. Here the overshoot is -194.22 arcsec, with the absolute error dropping below 
0.001 arcsec in 2.32 sec. 


• Results for Payload 2, angles about pitch axis (i.e.: “66") 


(arcsec) 


Absolute Pointing Error 1 
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Alternate Dissipative Controller 

The decentralized controller described previously was found to give satisfactory perfor- 
mance for the space station example. However, coupling between rigid and elastic motion was 
not severe, as the masses and inertias of the payloads were small compared to the station values, 
and the station represents a nearly rigid spacecraft. Decentralized control does not guarantee 
stability for the general case of a highly flexible LFSP with large inertial payload/LFSP coupling, 
as control system interactions can be destabilizing. An alternate method is to utilize a dissi- 
pative controller which requires that the matrix “A” in the second order vector equation (33) 
be symmetric. Development of the symmetric form of the equations of motion may be found in 
Appendix B. By redefining the system variables “brf 1 as in (35), it can be shown that the resulting 
“A” matrix is always positive definite. Using the sensor outputs as in (36) and the control law as 
in (37), it can be proven that the system is asymptotically stable (in the sense that all rotational 
and elastic motion tends to zero as time goes to infinity). The stability is guaranteed regardless 
of errors in the system parameters or ignored (unmodeled) elastic motion. Research is presently 
in progress on this type of controller. 


• Symmetric Form of Equations 

A6rj + B<5?) + C Srj = E£u 


where: 


A T = A > 0 
B T = B > 0 
C T = C > 0 


E = 


O 3 x3 


hx 3 

^0o(nx3) 

0[(6+n)x3n p ] 

03x3 


0.3x3 

I(3rip X3n p) 


8r]= (<5R£ 0 , <5aJ, *q T , 6a.J 0 , ... , <5a^ p0 ) 


y = E T r ? = 


[station attitude (rigid + flexible), at CMG location 
transformed gimbal angles, payload 1 


transformed gimbal angles, payload n p 

• Control Law 

6u = -G p 8y - G r 6y ; Gj = G p > 0 , Gj = Gj > 0 

• Closed-loop system asymptotically stable, i.e.: 

(£a<„ <5q, <5ajo, <5a 2o , ■ . ■ , <5a np0 ) -+ 0 as t -* «> 


• Research presently in progress 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 
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Concluding Remarks 


The problem of dynamic modeling and control of composite systems consisting of a central 
flexible space platform and articulated rigid payloads was considered. Generic linearized equations 
of motion were derived for a single rigid payload mounted on such a platform using a three-degree- 
of-freedom gimballing mechanism. This generic large flexible space platform (LFSP) model was 
then extended to the case of multiple independently pointed payloads. A simple decentralized 
control law was proposed for platform attitude control and payload pointing. Simulation results 
were obtained for an example problem with a single payload, and also with two payloads, attached 
to the space station. The results demonstrate the effects of dynamic couplings in the composite 
system, and also indicate that the control law used provides satisfactory payload pointing and 
platform stabilization for the example problem. 

An alternate centralized dissipative control law, which uses a symmetric form of the equa- 
tions of motion, was also proposed. This control law guarantees stability regardless of modeling 
erros and unmodeled modes. Further research is in progress on that topic. The model obtained 
herein is mathematically tractable, and yet has an accurate structure that includes all internal 
dynamic couplings. This model offers a suitable tool for analytical and numerical investigation 
of the dynamics and control of an important class of missions arising in the Control-Structure 
Interaction (CSI) program. 


• Developed nonlinear math models for LFSP/ articulated payload system 

• Extended to multi-payload systems 

• Obtained linearized equations 

• Simulation results for an example problem 

- indicates satisfactory performance for decentralized control law 

• Proposed a centralized dissipative control law 

- provides robust stability 

- development in progress 

• Further investigation: 

- extension to flexible articulated appendages 

- automated “FEM-type” model assembly 

- incorporation in CSI optimization problem 
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Appendix A: Equations for Multi-Payload Case 


The equations for an LFSP supporting multiple payloads will be derived from the 
linearized equations (8) through ( 12 ). Here a large flexible space platform described in part 
by n known elastic modes, is acted on by Jif external forces and tim moments exclusive of 
the torques required to point the n p payloads. 

The translational expression is: 

Tl p n p 

m 0 6R no = Y^6 F oi - Y ^ S Fi (-41) 

*=1 1 = 1 

The expression for the rigid-body rotation of the LFSP becomes: 

n m np n p n P 

ia = £ + E r t/ F oi - E - E {T i t- 42 ) 

1==1 t = 1 : 1=1 

Again the “prime” notation with RE and rE indicates the (3 x 3 ) matrix cross-product form 
of the vectors R 0 i = o^oi + t/,62 + 2,63 and r 0 j = ^oi^i + Voi ^2 + ^01^3 ■ The linearized 
equation for the elastic motion of the platform is: 

TIM 71 F U P W P 

6q + D<5q + A<5q = ^ ^ $ U oi£F oi - ^ (A 3 ) 

1 = 1 l—l 1=1 *=1 


Next are the rigid payload equations. There will be a translational and a rotational 
equation for each of the rip payloads. The CM of the payload (“point i* ) from the inertial 
reference for the deformed spacecraft is given by: 

R n i* = R no + Roi + + r i* (^) 

Differentiating this equation twice results in the inertial acceleration of the payload CM 
in the form of equation ( 7 ). Noting that the only external force acting on the payload is the 
constraint force, and linearizing gives: 


m ! {^R„o-Ki + (Cj i r i .)']« 5 a ( 


+ [* 5 -( c 2 r ,.)'* S ]^ q - 

-(Cj i r i .) , Cj i 6a i }=6F : 


(AS) 


C e ; is the nominal transformation between the elastic and some i th payload axes. The payload 
rotational equation of motion becomes: 


Iio(*So + + ClSaa) = 6T-, - (C5n.)'«F, 


(AS) 
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Here Ii 0 = C J- I; C e i is the i th payload inertia matrix transformed into the o-system. This 
inertia matrix multiplies the absolute inertial angular acceleration of the i th payload. 


Using summations over the n p payloads, equation (A5) can be used to remove 
terms from the LFSP equations (Al, A2, and A3). Furthermore, SF\ may be removed from 
the i th payload rotational equation (A6), yielding n p additional equations. This system of 
equations may then be assembled in the following matrix form: 

A Srj + B Si] + C Sr] = E<5u (A7) 

The vector Srj is dimensioned (6 + n + 3 n p ) x 1, and is given by: 

Sr/ — 8( X, Y, Z ; (f> 0} 9 oy ijj 0 ] gq, ... , q n \ <f )\ , $i, ^ 1 ; ... ; <^ np , 9 Upy xjj np (AS) 

The vector of force and moment inputs 6u is dimensioned (3 uf + 3tim + 3n p ) x 1, and is 
written as follows: 


<5u = 6( F o1 t , . . . , F onr T ; M 0 i T , . . . , M 


T, m T 

onjvi i ? 


T 

> -A-n, 


T nT 


(A9) 


where 6F 0 j, £M 0 ;, and 6T; are external LFSP forces and moments, and payload gimbal 
torques, respectively. 


The coefficient matrices A, B, and C, dimensioned (6 + n + 3 n p ) square, and E, 
dimensioned (6 + n + 3n p ) x (3 np + 3 njw + 3 n p ), are given below. Matrix A is essentially a 
“mass/inertia” matrix given by: 



Al[(6+n) x(6+n)] 

A2[(6+n)x3n p ] 

A = 

. A 3 [ 3n?)>< ( 6 + n )] 

A4(3n p X 3n p ) 


The blocks comprising matrix A are: 


Ax 



Tip 

+ Y m { 


1=1 


L 3X3 


n P 

- Y m,R! 


i=l 


Y i 
1= 1 


Y m > R oi - Y Y ^«Roi^i 


i=l 

i=l 

1=1 

n p 

Tip 

Tip 

Y m &i 

t= 1 

- Y 

X- 1 

Inxn + Y m 




-m np r' np Cj np 

a 2 = 

m iRoi ? i c Ji ■ 

Ron p ^rip C er)p 





(A10) 


(Alla) 


(All 6) 
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As — 


r f x mi Ii - mi^Ri Ii&ex + 


r' m. 


In p m n P Y n u ^n p 9n 0 ~b m n p r np ^n p J 


(411c) 


where: 


A 4 = 


(!i ~ m 1 r' 1 ri)Cj 1 


0 


3x3 


03XJ 


■" (I.. P - m », r n p < p ) C S 


FI = (Cjr,. )' R.; = (Roj + rj) 4>i = (#5 - 


(And) 


(A 12) 


Matrices B and C are for the most part null except for diagonal blocks multiplying 
the n modal coordinates. B and C are then modal damping and modal stiffness matrices 
written as: 

B = diag[0(5 x6) j D( nXn ), 0(3 npX 3 np )] 

C — diag[0( 6x6 ), A( nXn ), 0(3n p x3n p )] 


where: 

D = [diag(2(^)] (nxn) 
A = [diag(a;?)] {nxn) 
The coefficient matrix E is written as: 


(A13) 


(414) 


E = 


m) x 3rtjp] ^2[(6+n) x3tim ] ^3[(6-|-n) X3n p ] 


[_ 0(3n p x3nj7) f^(3n p x3nM) I(3n p x3n p ) 


(Alb) 


where the first row of blocks in E are given by: 



" I 3 X3 

Iax3 


’O 3 X 3 

0sx3 


"03x3 

03 x 3 

E x = 

r ol ’ 

. r ; 

' * * onp 

, E 2 = 

13x3 

I 3 X3 

, E 3 = 

— 13x3 

— 13x3 


- *uol 

<J> 

• ■ * uonp J 


-^<9ol 



- —&01 

- 


(A1S) 
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Appendix B: Symmetric-Form Equations for Single and Multi-Payload Cases 


A symmetric form of the composite system equations of motion will now be derived 
from the linearized asymmetric form. The linearized LFSP translation, rotation, and elastic 
equations of motion, as well as the payload rotational equations follow: 


(m 0 + m p )I 3 X 3 < 5 Rno - m p K'6a 0 + m p $6q - m p r'Cj p Sa p = ^<5F oi (B 1) 


1=1 


m p R4 p 6R no + [I 0 - m p Ro p R']6a 0 + m p R' op <f>5q - m p R; p r'C ep 6a p 

Tip ( JJO) 

= E r ^ F “+E <M ”>- iT p 

1=1 i = 1 


m p ® u 5R no - m^uR'^ao + [I nXn + m p $ u &]6q - m p $ u r'Cj p «5a p 


n f 


n\f 


^ $uoi<5F oi + # eo i^M oi - 


1=1 


1=1 


(B 3) 


m p r'<5R no + [I po - m p r'R']<5a 0 + + m p r'#](!q 


+ (I po - m p f'r')Cj p ia p = <5T p 


(B 4) 


where: 


-/ _ /pT „ y 1 is/ _ rryf » ^ _ ( x * T 

1 — V^ep r pC 


R’ = (R' +r') * = (*I-r'*J) 


(B 5) 


A symmetric form of the equations of motion may now be obtained by removing <5T p 
from equations (B2) and (B3) using (B4), yielding: 


ra p R'£R no + (I 0 + Ipo - m p R'R')6a 0 + (I po 3»J + ?n p R'$)6q 

+ (Ipo - m p R'r')Cj p 6a p = £ r^F.,, + £ «M t 


jjm 


(Be) 


i= 1 


i= 1 
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(B7) 


m p ^ T <5R no + (^^Ipo — mp$ T R/)£a 0 + (I nX n + 3»0l po ^J + m p <l» T 3>)<5q 

™F tim 

+ ($ e Ipo ~ m p ^ T f'C^ p 6a p = ^ $uoi<5F 0 i + $<? 0 ;6M oi 

i=l »=1 


A matrix form of equations (Bl, B 6 , B7, and B4) may be assembled as follows: 


A Srj + ~B6rj + C 8rj — E^u (£ 8 ) 

Here A is a (4 x 4) block matrix of overall dimension (n + 9) square. Matrices B and C are 
essentially null except for n diagonal terms corresponding to the elastic motion. The vector 
8i 7 of system parameters is (n + 9) x 1, and is given by: 


6r ] = 6( RL;aJ;q T ;aJ 0 ) 


(B 9) 


The transformed payload angles £a po = Cj p £a p must appear in the vector of system param- 
eters for symmetric equations to result. The vector of control inputs is: 


iu = («Fj 1 , 


..,«F 




T\T 


(BIO) 


The equations for an LFSP supporting multiple payloads will be given as obtained 
from the linearized equations (Bl, B 6 , B7, and B4). Here a large flexible space platform 
described in part by n known elastic modes, is acted on by rip external forces and um moments 
exclusive of the torques required to point the n p payloads. 

The equations of motion for multi-payload configurations may be written in the fol- 
lowing form: 

ASrj + B<5?7 + CSrj = E<5u (£11) 

The vector 8rj is dimensioned (6 + n + 3 n p ) x 1 , and is given by: 

6n = 6( RL;a?;q T ;aL ... ,a^ p0 ) T (£12) 

The vector of force and moment inputs 6 u is dimensioned (3 np + 3n^ + 3 n p ) x 1 , and is 
written as follows: 

8u = 8(F o1 , ... , F onF ; M o1 , ... , M onM ; T l5 ... , T„ p ) (£13) 

where £F 0 ; , <5M 0 ;, and <5T; are external LFSP forces and moments, and payload gimbal 
torques, respectively. 
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The coefficient matrices A, B, and C, dimensioned (6 + n + 3 n p ) square, and E, 
dimensioned (6 + n + 3 n p ) x (3 n F + 3 um + 3 n p ), are given below. Matrix A is essentially a 
“mass/inertia” matrix given by: 


A = 


-^•l[(6 + n) x(6+n)] -^-2[(6+n) x3n p 

L •^■3[3n p X (6+n)] A 4 ( 3npX3nj) ) J 
The blocks comprising matrix A are: 


(314) 


Ai 


m c 




r \ — 

o + E m i] : 3X3 - E m i R [ 


E 


E m i R i 


i= 1 


E ">,*? 
1=1 


i = 1 


i= 1 


Io - E (Iio - m< R | R 0 E (!io + m,R!$i) 
1=1 1=1 

- E (««Ilo - R -0 Inxn + E 


(315a) 


i=l 


+ E">i*?*i 

i= 1 


A o = 


-mir'i 


I lo - mj R^r^ 


“ m n P r n P 


In p o m n p R n p r n p 
^fln p ln p o ^*ip^n p **n p 


(3156) 


mifi 


Iio - m^Ri. Iio^ + miri*! 


m n p r n p In p o m n p r n p R n p ln p o^ 0 np 4” m «p r n p ^ n P -J 


(315c) 


where: 


A 4 = 


Iio - mir^ri ... 


0 


3x3 


0sx3 


... I n o — m n I*' r' 

n P° 7l p n p n p 


rj = (Cj,r,. )' ft! =(Ri, + r|) Si = (*S-^*5) 


(B15d) 


(B 16 ) 
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Matrices B and C are for the most part null except for diagonal blocks multiplying 
the n modal coordinates. B and C are then modal damping and modal stiffness matrices 
written as: 


B diag[0( 6x6 ), D( nXn ) ? 0(3u p x3n p )] 

C diag[0( 6x6 ), A( nXn ), 0^3 npX 3 np j] 


(B 17) 


where: 

D = [diag(2£u;j)](„xn) 
A = [diag(w, 2 )](nxn) 


(BIS) 


The coefficient matrix E is written as: 


E = 


E l[(6-f-n)x3n F ] 
®(3n p x3nip) 


E 2[(6+n)x3n M ] 
®(3 n p x3«m) 


0[(6 + n)x3n p ] 
I(3n p x3n p ) 


where Ei and E 2 are given by: 



" ^3x3 

I 3 X 3 


"03x3 

03x3 

Ei = 


r 7 

onp 

II 

w 

I 3 X3 

13x3 


- $uol 

tf> 

■ • * uonp J 


-**ol -■ 

■ • ^ ^onjvj - 


(519) 


(520) 
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